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ABSTRACT 


The  simulation  of  a  plasma  sheath  using  an  artificial  dielectric 
is  studied  in  this  investigation  and  is  applied  to  an  antenna  geometry 
which  is  similar  to  some  configurations  encountered  in  aerospace 
applications.  The  antenna  configuration  is  equivalent  to  a  horn  in 
an  infinite  ground  plane  with  an  unbounded  plasma  layer  in  front  of 
the  horn.  The  plasma  layer  is  simulated  by  a  rodded  medium,  and 
the  radiation  patterns  of  the  antenna  system  are  studied  experimentally 
at  y,  10,  and  11  KMC. 

The  characteristics  of  a  rodded  medium  are  obtained  both  theo¬ 
retically  and  experimentally.  It  is  shown  that  a  rodded  medium  not 
only  simulates  the  isotropic  characteristics  of  a  plasma  sheath  but 
also  the  anisotropic  characteristics  of  a  plasma  in  an  infinite 
steady  mtgnetic  field.  If  the  rods  are  oriented  in  one  directio  , 
the  equivalent  tensor  permittivity  of  the  rodded  medium  is  in  the  same 
form  as  the  tensor  permittivity  of  a  cold  plasma  in  a  very  strong 
magnetostatic  field.  The  theoretical  analysis  is  based  on  the 
geometrical  symmetries  of  the  medium  and  the  transmission  line 
approach  to  propagation  in  an  artificial  dielectric  of  periodically 
spaced  conducting  obstacles. 

The  anisotropic  characteristics  of  a  rodded  medium  are 
demonstrated  in  X-band  waveguide  measurements. 

Two  constructions  of  rodded  media  are  utilized  in  the  pattern 
measurements.  The  experimentally  determined  radiation  patterns 
seem  to  have  u  secondary  maximum  associated  with  the  properties  of  the 
plasma  sheath.  More  energy  is  radiated  in  the  off-axis  region  when 
the  plasma  sheath  is  present. 
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I.  INTRODUCTION 


The  purpose  of  this  investigation  was  to  study  plasma  simulation 

using  a  rodded  medium  and  to  study  the  radiation  patterns  of  an  X-band 

horn  antenna  in  the  presence  of  a  simulated  plasma  sheath.  The  plasma 

sheath  is  assumed  to  be  under  dense  or  operating  above  the  plasma 

frequency.  The  antenna  geometry  is  shown  in  Figure  l.  Several 
12  3  4 

authors  *  '  '  have  studied  this  geometry  or  very  similar  geometries 
previously,  but  to  obtain  a  clean  experiment  using  a  laboratory  plasma 
is  very  difficult.  Typically,  the  experimental  and  theoretical 
geometries  are  not  the  same.  Usually  the  theoretical  models  arc 
highly  idealized  and  only  approximate  the  actual  experimental  geometry 
because  of  the  inherent  difficulties  in  obtaining  a  solution  to 
Maxwell's  Lquations  for  the  actual  geometry  of  interest.  For  example, 
the  theoretical  models  pertaining  to  this  study  assume  that  the  plasma 
sheath  is  unbounded }and  the  ground  plane  is  infinite  in  extent,  while 
the  physical  dimensions  of  the  laboratory  apparatus  are  finite  in 
size.  The  glass  walls  which  confine  the  gaseous  discharge  may  also 
complicate  the  experimental  geometry.  The  theoretical  models  neglect 
nonuniformities  in  the  sheath  which  may  be  present  in  the  actual 
laboratory  experiment.  In  the  final  analysis  the  experimental 
measurements  prove  or  disprove  the  validity  of  the  theoretical  model 
used  to  approximate  the  actual  geometry  under  investigation. 

Theoretical  Models 

7 

Umura1"  obtained  the  expressions  governing  the  far  field  radiation 
patterns  of  an  infinitely  long  slot  antenna  in  a  ground  plane  which  is 
covered  with  a  plasma  sheath.  The  analysis  is  based  on  Fourier  Trans¬ 
form  techniques  and  assumes  two  different  aperture  distributions. 

In  the  evaluation  of  the  far  field  patterns  the  contributions  of  surface 
and  complex  waves  are  neglected.  One  ot  the  distributions  is  equal  to 
a  constant  electric  field  developed  across  the  gap.  The  other  distri¬ 
bution  assumes  that  the  electric  field  is  parallel  to  the  edge  of  the 
slit  and  has  a  co-sinusoidal  distribution  across  the  gap.  The  b  plane 
and  II  plane  radiation  patterns  of  the  horn  geometry  can  be  approximated 
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APERTURE 


a  =3.25  inches 
b  =2.625  inches 
1=  3  to  S.Oinches 


Antenna  geoaetry  used  in  experiaentul  study 


by  the  two  equivalent  slit  problems0  described  above.  Aside  from  the 
cylindrical  or  spherical  wave  terms  of  the  far  field  radiation  patterns, 
the  algebraic  expressions  which  relate  the  relative  pattern  changes 
of  the  horn  and  equivalent  slit  in  the  presence  of  the  plasma  are 
approximately  equal  and  identical  when  the  plasma  layer  is  removed, 
ihe  radiation  from  the  slit,  which  has  a  constant  electric  field 
developed  across  the  gap,  is  given  by  Equation  (1)  and  is  equal 
to  the  E  plane  pattern  of  the  horn  geometry. 
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Similarly,  the  patterns  obtained  from  the  slit,  which  are  excited  with 
the  co-sinusoidal  distribution,  resemble  the  H  plane  patterns  of  the 
horn  antenna  and  arc  expressed  by  Equation  (2) . 


2 

n 

T“ 


-t2  r 


1 


- 5 — 

2/-..  „•»  /  cos  ? 

c°s  (11)  ♦  -J - -j  , - , 

n  -  sm  6 


sin^K  i) 


J(2) 


The  index  of  refraction  of  the  plasma  is  given  by 

2 

n2  -  1  -  , 

at 

where  is  the  plasma  frequency  and  w  is  the  wave  frequency.  The 
dimension  of  the  aperture  and  the  thickness  of  the  plasma  layer 
are  equal  to  a,  b,  and  £.  respectively  .The  equivalent  Z  component 
of  the  wavenumber  in  the  plasma  is  equal  to 
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Flock  *  studied  the  horn  geometry  from  both  a  theoretical  and 
experimental  standpoint.  The  theoretical  model  assumes  that  the 
dominant  pattern  changes  are  caused  by  the  currents  induced  in  the 
plasma  region  directly  in  front  of  the  horn  and  all  other  currents  in 
the  plasma  are  zero.  The  currents  directly  in  front  of  the  aperture 
are  assumed  to  have  a  co-sinusoidal  distribution  in  the  X  direction 
similar  to  the  distribution  of  magnetic  current  in  the  aperture. 

The  analysis  neglects  the  effects  produced  by  the  currents  which  are 
induced  in  other  regions  of  the  plasma  by  the  near  fields  of  the  horn. 
The  theoretical  calculations  predict  that  the  radiation  from  the 
equivalent  source  in  the  plasma  will  have  a  maximum  at  an  angle  6^ 
given  by 


cos  0.,  *  n  , 

M 

when  u  is  greater  than  u>  . 

P 

The  gaseous  discharge  which  he  used  in  the  experimental  apparatus 
had  dimensions  of  the  order  of  an  X-band  standard  gain  horn  which 
restricts  the  interpretation  of  the  experimental  results  in  terms 
of  the  unbounded  plasma  layer  model. 

Artificial  Dielectrics 

The  radiation  patterns  described  in  this  report  have  been  studied 
experimentally  using  an  artificial  dielectric  to  simulate  the 
properties  of  the  plasma  sheath.  An  artificial  dielectric  is  used 
to  alleviate  some  of  the  experimental  difficulties  of  performing 
the  measurements  with  a  laboratory  plasma.  The  artificial  dielectric 
slab  can  be  made  large  enough  to  approximate  the  conditions  of  an 
unbounded  plasma  layer.  Also,  the  electromagnetic  properties  of 
an  artificial  dielectric  are  homogeneous  and  are  characterized  by  an 
equivalent  plasma  frequency  which  can  be  determined  theoretically 
within  reasonable  limits.  The  experimental  measurements  can  be 
performed  in  a  steady  state  condition  and  do  not  have  to  rely  on 
the  shot  to  shot  reproducibility  of  a  gaseous  discharge.  The  main 
attribute  of  the  artificial  dielectric  is  the  ease  of  controlling 
the  experimental  variables  such  as  the  size  of  the  slab,  the 
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equivalent  plasma  frequency  or  the  thickness  of  the  sheath. 

Recently,  it  was  observed**  that  the  complex  dielectric  constant  of 
a  homogeneous  isotropic  plasma  could  be  simulated  by  a  rodded  medium. 

The  currents  induced  in  the  rods  tend  to  simulate  the  oscillating 
inductive  currents  in  a  plasma  thereby  making  the  simulation  possible. 

In  general,  the  rods  form  a  cubical  lattice  if  the  dielectric  properties 
are  to  be  isotropic  for  all  polarizations. 

Like  a  crystal  lattice  structure,  an  artificial  dielectric  can 

7  8  9 

also  exhibit  anisotropic  '  *  behavior  under  certain  conditions.  In 
principle,  it  should  therefore  be  possible  to  simulate  the  anisotropic 
behavior  of  a  plasma^  in  an  infinite  magnetostatic  field  with  a  rodded 
medium,  for  these  conditions  the  electrons  in  the  plasma  are  confined 
to  move  along  the  magnetic  lines  of  force.  Anisotropic  characteristics 
of  a  rod  medium  are  observed  when  the  rods  are  oriented  in  one  or  two 
directions.  The  induced  currents  are  restricted  to  flow  along  certain 
preferred  directions,  namely  along  the  rods.  Therefore,  the  propagation 
characteristics  of  the  rodded  medium  are  a  function  of  the  polarization 
of  the  electric  field  making  the  medium  anisotropic. 

The  technique  used  in  this  report  can  be  useful  in  the  design  and 
evaluation  of  re-entry  vehicle  antenna  systems.  The  rodded  medium  can 
be  scaled  to  simulate  an  over  dense  as  well  as  an  under  dense  plasma 
sheath.  An  artificial  dielectric  can  not  simulate  all  the  complexities 
of  a  re-entry  plasma  sheath,  but  at  this  stage  the  simulation  of  the 
simplest  properties  of  a  plasma  can  be  very  useful  in  laboratory  and 
ground  testing.  A  great  deal  of  insight  into  the  problems  of  re-ent^y 
communications  can  be  gained  from  these  simplified  model  studies,  and 
it  is  advantageous  to  understand  the  behavior  of  the  simplified 
geometries  before  attempting  to  solve  the  more  realistic  problems, 
idany  times  the  simplified  analysis  can  be  used  as  a  guide  in  the  more 
complex  geometries. 
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Dielectric  Constant  of  a  Plasma 

The  electromagnetic  properties  of  a  plasma  can  be  represented 

11  12 

by  a  complex  dielectric  constant  which  has  tensor  properties 

under  the  application  of  a  static  magnetic  field.  For  the  purpose 
of  this  discussion,  the  plasma  is  assumed  to  be  a  homogeneous 
slightly  ionized  gas  which  is  electrically  eutral.  Collisional 
effects  are  assumed  to  be  zero  and  the  electrons  haye  no  thermal 
motion.  The  dielectric  constant  for  the  plasma  is  determined  by 
solving  Maxwell's  Equations  in  conjunction  with  the  equation  of 


motion  of  the  electrons. 

V  x  E  =  -  jw  y  11  (3) 

o 

V  x  11  *  jw  e  E  ♦  J  (4) 

J  *  p  v  (5) 

m  »  e  [E  ♦  v  x  B]  (6) 


The  electron  charge  density  is  equal  to  p  and  the  velocity  of  the 
electrons  is  equal  to  v.  The  curl  of  K  can  also  be  expressed  in 
terms  of  a  tensor  permittivity. 

V  x  li  *  jw  [c]  [E]  (7) 

If  the  above  equations  are  linearized,  the  resulting  tensor 
permittivity  is  given  by 

ell  je12 
-jei2  e22 
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when  th$ . .  magnetic  field,  B  ,  is  applied  in  the  X,  direction. 

w 

The  rectangular  components  of  the  tensor  are  equal  to 


B  e 

wc  ■  -  -jjj—  »  electron  cyclotron  frequency 


if  collisional  effects  arc  neglected.  It  should  be  noted  that  e 
is  a  negative  quantity  for  an  electron  and  is  equal  to  the  electronic 
charge.  N  is  the  electron  density,  tQ  is  the  permittivity  of  free 
space,  and  m  is  the  mass  of  an  electron. 


If  the  static  magnetic  field  becomes  infinite  in  strength,  the 
tensor  [ef]  simplifies  to 
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Similarly,  when  Bq  goes  to  zero,  the  plasma  becomes  isotropic  with 
a  dielectric  constant  given  by 


er  " 


w  * 

i  -  4-  ; 
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B 
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0 


It  will  be  shown  that  the  above  behavior  is  simulated  by  a 
rodded  medium  when  there  are  many  rods  per  wavelength. 
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II.  RODDED  MEDIUM 


The  dielectric  properties  of  the  roddcd  medium  will  be  studied 

in  this  section.  The  rods  are  assumed  to  be  perfect  conducting 

cylinders  therefore  making  the  equivalent  collision  frequency  zero. 

The  electromagnetic  properties  of  the  rodded  medium  can  be  described 

by  a  macroscopic  dielectric  constant  which  has  tensor  characteristics. 

If  the  rods  are  only  oriented  in  one  direction,  the  roc'ded  medium 

exhibits  anisotropic  behavior, and  the  propagation  constant  depends 

on  the  polarization  of  the  electric  field.  If,  however,  the  rods 

form  a  cubical  lattice,  the  dielectric  properties  become  isotropic. 

The  spacing  between  adjacent  rods  must  be  much  smaller  than  a  free 

space  wavelength.  As  the  wavelength  approaches  the  rod  spacing, 

diffraction  effects  become  dominant, and  the  rodded  medium  no  longer 

acts  like  a  dielectric  material.  Also,  the  air  artificial  dielectric 

13 

interface  behaves  incorrectly  when  the  rod  spacing  becomes  an 
appreciable  portion  of  a  wavelength. 

The  equivalent  tensor  dielectric  constant  for  an  artificial 
dielectric  with  rods  in  the  X^  direction  only  will  be  determined. 

It  can  be  shown  from  symmetry  arguments  that  the  tensor  permittivity 
of  the  rodded  medium  has  the  following  form  (see  Appendix  B). 


[c]  *  e 
1  1  o 


The  rodded  medium  can  be  represented  by  an  infinite  number  of 
periodically  spaced  wire  gratings  when  propagation  along  the 
or  X^  axis  is  considered.  This  analogy  leads  to  the  equivalent 
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transmission  line  circuit^*  15,16  £or  j^ded  oedimn  illustrated 
in  Figure  2.  The  surface  impedance  Z  represents  the  energy 
storage  in  the  vicinity  of  the  wire  gratings.  If  the  electric 
field  is  polarized  parallel  to  the  rods,  the  energy  near  the 
gratings  is  stored  in  the  magnetic  field  making  the  equivalent 
surface  impedance  inductive.  When  the  electric  field  is  perpendicular 
to  the  wires,  the  equivalent  surface  impedance  is  capacitive.  The 
resulting  propagation  constant  of  the  medium  therefore  depends  on 
the  polarization  of  the  electric  field. 

The  equation  governing  the  propagation  constant  of  a  periodically 

17  18 

loaded  transmission  line  *  with  shunt  impedances  is  given  by 

Z 

cos  Kd  »  cos  kQd  ♦  j  sin  k  d  ,  (14) 

s 

where  K  is  the  average  propagation  constant  of  the  rodded  medium, 

k^  is  the  free  space  wave  number,  ZQ  is  the  characteristic  impedance 

of  free  space,  and  Z  is  the  equivalent  surface  impedance  of  the 

wire  gratings.  First,  assume  the  electric  field  is  polarized  in 

19 

the  direction.  The  surface  impedance  Z$  is  equal  to 


where  d  is  the  distance  between  rods,  r  is  the  radius  of  the  rods 
and  \Q  is  the  free  space  wavelength.  The  subscript  | j  denotes  that 
the  electric  field  is  parallel  with  the  wires.  Equation  (15)  is 
valid  when  d  is  less  than  Aq/2,  and  r  is  much  less  than  d.  When  <£  is 
equal  to  or  greater  than  Xq/2,  the  plane  of  wires  starts  to  take 
on  the  character  of  an  optical  grating  and  many  higher  order 
reflected  beams  are  sometimes  observed.  Propagation  characteristics 
along  either  the  Xj  or  X^  axis  are  therefore  determined  by 


-10- 


EQUIVALENT  GRATING  STRUCTURE 
OR  SURFACE  IMPEDANCE 


EQUIVALENT  TRANSMISSION  LINE 
ANALOGY  FOR  PROPAGATION 
ALONG  THE  X,  OR  X2  AXIS. 


Redded  nediua  geometry 


< 


I 


Equations  (14)  and  (15) 


cos  K  d 

II 


cos  k  d 
o 


sin  k  d 
o 


(16) 


so  long  as  the  electric  field  is  parallel  to  the  axis.  Equation 
(10)  can  be  simplified  if  kQd  is  assumed  to  be  much  less  than 
unity.  Then  the  above  expression  deduces  to 


1 


2*  d  In 


(17) 


The  relationships  between  the  propagation  constant  K 


II 


tensor  dielectric  constant  [e]  are  shown  in  Table  (la) 

2 

For  the  above  conditions,  the  propagation  constant,  Kjj 
equal  to 


and  the 
in  Appendix  A. 
,  is 


(18) 


The  combination  of  Equations  (17)  and  (18)  yields 


e33 


(19) 


If,  however,  kQd  is  greater  than  0,1,  should  be  determined 
from  Equations  (16)  and  (18). 

Now,  assume  that  the  electric  field  is  parallel  to  axis 

and  propagation  is  in  the  X  direction.  For  this  polarisation 

^  20 

the  equivalent  surfaco  impedance  is  given  by 


\ 


12- 


(20) 


The  subscript  X  denotes  that  the  electric  field  is  perpendicular 
to  the  wires.  Lquation  (20)  is  valid  when  the  ratio  r/d  is  much 
less  than  unity.  The  equation  governing  the  propagation  constant 
for  tills  polarization  is  obtained  by  substituting  Lquation  (20) 
into  Lquation  (14). 


cos  Kid  *  cos  k  d 
o 


(2ht) 

mr 

o 


sin  k  d 
o 


Lquation  (21)  reduces  to 


,  2ttT 
l*  “ 


when  k^d  is  much  less  than  unity.  The  value  of  KjJ*'  is  again  obtained 
from  Tabic  (la)  in  Appendix  A 


and  is  combined  with  Lquation  (22)  to  yield 


It  is  interesting  to  note  that  Lquation  (24)  can  be  deduced 

2 1 

by  using  the  concept  of  electric  polarizability  ot  the  rods*  .  This 
approach  is  equivalent  to  the  techniques  used  in  determining  the 
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22  2  3 

dielectric  constant  of  a  molecular  medium  *  . 

For  most  practical  cases,  the  ratio  r/u  is  extremely  small, 

-2  -3 

typically,  ranging  from  10  to  10  .  Under  these  conditions, 

liquation  (24)  reduces  to 

£„  %  1  ;  r2/d2  «  1  ,  (25) 


In  a  similar  manner,  e,,-,  is  also  equal  to  unity.  The  resulting 
dielectric  tensor  for  an  artificial  dielectric  with  rods  in  the 
Xj  direction  is 


[c] 


1 

0 

0 


1 

0 


0 


0 


2  2 

1  * 

o  p 


(20) 


where 


(27) 


When  the  ratio  d/XQ  is  between  0.02  and  0.5,  Equations  (16)  and 
(18)  can  be  used  to  determine 

Equation  (26)  resembles  the  tensor  dielectric  constant  of 
an  anisotropic  plasma  with  an  infinite  static  magnetic  field 
applied  in  the  direction. 


1 


0 


0 


[c]  = 


2.2 

i  -  ui  /u 

P 


(28) 


The  of  the  roddcd  medium  is  therefore  the  analogue  of  the 
cutoff  wavelength  in  the  plasma. 

When  the  plasma  is  subjected  to  an  infinite  magnetostatic 
field,  the  oscillating  electron  currents  in  the  plasma  are  con¬ 
fined  to  flow  along  the  magnetic  lines  of  force.  In  a  similar 
manner,  the  oscillatory  currents  in  the  rodded  medium  are  con¬ 
fined  to  flow  along  the  wires  thereby  simulating  the  effect  of 
applying  an  infinite  steady  magnetic  field  to  the  plasma. 


The  tensor  dielectric  constant  for  an  artificial  dielectric 
having  rods  parallel  to  both  the  and  X^  axis  is  given  by 

r  i 


\  i 

i  * 


0  1  -  \2/x2 

o'  p 


1  *  Xo/Xp2 


(29) 


Similarly,  the  tensor  dielectric  constant  for  a  cubical  array 
of  rods  is  equal  to 


d/d 

o  p 


0  1  - 


.2.2 

A  /X 

o  P 


!  2  2 
1  -  A  /A 
o  p 


(30) 
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Since  all  the  diagonal  terms  are  equal,  the  medium  is  isotropic  and 
the  dielectric  constant  is  given  by 


c 


-x2/x2 

o'  p 


(31) 


The  equivalent  index  of  refraction,  n,  for  an  isotropic  rodded 
medium  is  obtained  from  an  alternate  form  of  bquation  (16). 


cos 


(32) 


The  index  of  refraction  of  the  medium  is  defined  by 


n  *  Ao/A  , 


(33) 


where  X  is  the  average  wavelength  in  the  rodded  medium.  The  index 
of  refraction,  as  given  in  Equation  (32),  can  be  approximated  by 


n 


1  - 


xl,xl 

o  p 


1/2 

* 


(34) 


where 

P 


is  the  largest  zero  of  F(A) 


F(A)  =  1 


0. 


(35) 


Equation  (35)  is  the  result  of  setting  n  equal  to  zero  in  Equation  (32). 

X  is  therefore  the  longest  cutoff  wavelength  of  the  rodded  medium. 

P 

When  the  spacing  constant  d  is  much  smaller  than  a  wavelength,  the 
equivalent  plasma  wavelength  reduces  to  the  value  given  in 
Equation  (27) . 
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,n  (rvr 


(36) 


X 


2 

P 


The  theoretical  values  of  n  for  the  rodded  medium  used  in  the 
experimental  study  are  shown  in  Figure  3.  The  pertinent  quantities 
pertaining  to  the  experimental  geometry  are  given  below 


d 

r 

freq  range 


0.3  in. 

0.002  in. 

8  to  12  kmc. 


In  general ,  the  rodded  medium  behaves  like  a  homogeneous  isotropic 
plasma  when  the  rods  form  a  cubical  lattice  structure.  In  certain 
idealized  geometries,  however,  the  cubical  lattice  structure  is  not 
necessary.  For  example,  the  reflection  and  transmission  characteristics 
of  a  plasma  slab  at  normal  incidence  can  be  simulated  by  orienting  the 
rods  in  the  direction  of  the  incident  electric  field.  Physically, 
rods  tend  to  simulate  the  effects  of  the  oscillatory  motions  of 
the  free  electrons  in  the  plasma.  If  any  of  the  rectangular  components 
of  the  electric  field  are  zero  throughout  the  plasma  region,  the 
rods  oriented  in  those  directions  can  be  eliminated  since  they  are 
not  involved  in  the  simulation. 
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FIGURE  3  -  Characteristics  of  the  rodded  Medium 


III.  EXPERIMENTAL  STUDY 


Waveguide  Measurements 

The  propagation  characteristics  of  a  rodded  medium  similar  to 

the  ones  used  in  the  pattern  measurements  have  been  measured  in  an 

X- band  waveguide  interferometer  shown  in  figure  4.  The  dominant 

TE^0  mode  in  a  rectangular  waveguide  is  analogous  to  two  uniform 

plane  waves  bouncing  off  the  walls  of  the  waveguide,  as  shown  in 

Figure  5.  The  equivalent  angle  of  incidence,  0,  of  these  waves  is 

24 

related  to  the  wav  guide  wavelength  and  is  given  by 

V*g  “  cos  0  •  (37) 


where  X  is  the  waveguide  wavelength.  The  equivalent  angle  of 

£  25 

incidence  can  also  be  given  in  terms  of  the  cutoff  wavelength 

of  the  X-hand  waveguide, 


X  /A  =  sin  6  , 
o  c 


(38) 


The  waveguide  measurements  are  equivalent  to  the  same  measurements 
performed  in  free  space  at  an  angle  of  incidence  given  by  Equation 
(38). 


The  relative  dielectric  constant  of  the  rodded  medium  can  be 
related  to  the  absolute  fringe  shifts  measured  with  the  inter* 
ferometer.  The  relative  fringe  shift,  Aarcj»  is  given  by 


ha 


rel 


NL 


-  e, 


(39) 


where  G^,  is  the  reading  on  the 
absence  of  the  sample  and  6^  is 
in  the  presence  of  the  sample. 


calibrated  phase  shifter  in  the 
the  reading  on  the  phase  shifter 
The  absolute  value  of  the  fringe 
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FIGURE  4 


Diagram  of  the  X«band  waveguide  interferometer 


equivalent  rodded  medium  interface 


CO 
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FIGURE  5  -  Rodded  medium  in  waveguide  geometry. 


shift  is  given  by 


Aa  *  Aarcl  1  M2ti  ’  (40) 

where  M  is  an  integer.  The  ambiguity  in  the  measured  value  of 
Aa  stems  from  the  inability  of  the  relative  phase  shift  measure* 
ments  to  distinguish  multiples  of  2«,  The  value  of  M  can  be 
determined  by  estimating  the  relative  dielectric  constant,  c  ,  or 
by  performing  an  additional  set  of  measurements  with  a  different 
length  sample. 

The  measured  value  of  Aa  can  also  be  related  to 

Aot  *  ^610  "  *  (41) 

where  0^  is  the  propagation  constant  of  the  empty  waveguide,  8Z 
is  the  propagation  constant  of  the  loaded  waveguide,  and  £  is 
the  length  of  the  sample.  Equation  (41)  is  only  approximate  in  that 
the  phase  shifts  attributed  to  the  multiple  reflections  off  the 
ends  of  the  sample  are  neglected.  When  the  magnitude  of  the 
reflection  coefficient  is  small,  the  measured  transmission  phase 
given  by  Equation  (40)  approaches  the  value  given  by  Equation  (41) . 

The  propagation  constant  of  the  loaded  waveguide  is  equal 


3 


z 


;  TEj0  mode 


(42) 


where  is  the  relative  dielectric  constant  of  the  sample,  and 

X  is  the  cutoff  wavelength  of  the  unloaded  waveguide.  Equation 
C 

(42)  cun  be  rewritten  in  the  form 


£ 


r 


♦ 


(43) 
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It  is  interesting  to  note  that  liquation  (43)  can  be  obtained  from 
Snell’s  Law  when  the  angle  of  incidence  is  given  by  Equation  (38), 


> 


which  again  illustrates  the  analogy  between  the  waveguide  and  the 
free  space  measurements. 

♦ 

The  results  of  the  measurements  performed  on  the  rodded  medium 
used  in  the  antenna  pattern  measurements  are  shown  in  Figure  3. 

The  measured  index  of  refraction  corresponds  to  the  case  with  the 
electric  field  parallel  to  the  wires. 

Additional  measurements  were  performed  with  the  wires  oriented 
perpendicular  to  the  electric  field.  In  one  case  -he  wires  were 
oriented  in  the  transverse  plane  but  perpendicular  to  the  electric 
I  field,  and  in  the  other  the  wires  were  oriented  in  the  direction 

of  propagation.  These  measurements  indicated  no  apparent  phase  shift 
due  to  the  presence  of  the  wires,  making  the  equivalent  index  of 
refraction  for  these  polarizations  equal  to  unity.  The  behavior 
of  a  rodded  medium  with  the  wires  oriented  in  one  direction,  there¬ 
fore,  exhibited  the  anisotropic  properties  described  in  the  previous 
section  when  viewed  in  the  X-band  waveguide  geometry. 

The  surface  impedance  of  a  single  plane  of  wires  was  measured 
with  the  aid  of  the  slotted  line  shown  in  Figure  4.  The  wires  were 
oriented  in  the  direction  of  the  electric  field,  and  the  experi¬ 
mental  data  and  theoretical  model  were  within  10%,  as  shown  in 
Figure  6.  Better  agreement  was  obtained  by  using  a  one  point  fit 
between  theory  and  experiment.  In  other  words,  the  equivalent 
surface  inductance  of  the  plane  of  wires  was  obtained  experimentally. 
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NORMALIZED  SURFACE  IMPEDANCE 
OF  A  PLANE  OF  WIRES. 
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^IGUR£  6  -  Surface  inpedance  of  a  single  plane  of  equally  spaced  wires. 


Antenna  Geometry 


The  radiation  patterns  from  an  X-band  horn,  terminated  in  an 
infinite  ground  plane,  which  was  radiating  through  an  unbounded 
homogeneous  plasma  sheath,  were  measured  experimentally  (see  Figure  1). 
The  sheath  was  assumed  to  be  under-dense  or  operating  above  the  plasma 
frequency.  The  actual  experimental  antenna  geometry  only  approxi¬ 
mated  the  idealized  model  described  above.  The  sheath  was  simulated 
by  an  artificial  dielectric  slab  which  was  large  compared  with 
the  dimensions  of  the  antenna  aperture  and  slab  thickness.  The 
restrictions  on  the  size  of  the  slab  were  a  necessary  requirement 
for  the  slab  to  approximate  an  unbounded  layer.  For  this  particular 
experimental  study,  the  aperture  size  was  given  by 

a  -  3.25  inches 

b  *  2 .625  inches, 

and  the  dimensions  of  the  ground  plane  were  about  six  times  greater 
than  the  dimensions  of  the  aperture.  The  axial  length  of  the  horn 
was  12  inches. 

The  horn  was  excited  with  a  TEjQ  mode  which  is  equivalent  to 
having  the  following  aperture  distribution 

E  (z  *  0)  *  u  E  cos  — ^ 

y  o  ft 

developed  in  the  mouth  of  the  horn. 

Design  of  the  Rodded  Medium 

Two  artificial  dielectric  geometries  were  utilized  in  the 
pattern  measurements.  Both  models  were  constructed  with  38  gauge 
wire  and  a  wire  spacing  of  0.3  inch.  The  equivalent  plasma 
frequency  of  the  artificial  dielectrics  was  8.1  KMC.. 


-25- 


The  majority  of  the  measurements  were  performed  with  the  rodded 

medium  shown  in  Figures  7,  8  and  9.  The  equivalent  thickness 

« 

of  the  slab  could  be  varied  from  0  to  3  inches  in  0.3  inch  steps. 

In  this  geometry  the  wires  were  oriented  in  the  Y  direction  and 
were  held  in  position  by  the  grooved  posts  located  at;  the  top  and 
bottom  of  the  ground  plane.  Wire  sag  was  eliminated  by  keeping 
the  wires  under  tension.  The  springs  located  near  the  upper  edge 
of  the  ground  plane  maintained  the  desired  tension  in  the  wires. 
Microwave  absorbing  material  was  placed  around  the  posts  and 
associated  hardware  to  eliminate  scattering  problems.  The  absorbing 
material  altered  the  measured  E  plane  patterns  at  angles  greater 
than  50  degrees,  but  these  effects  were  far  less  severe  than  the  ones 
occurring  when  the  absorbing  material  was  absent.  The  H  plane 
patterns  were  practically  unaltered  by  the  presence  of  the  absorber. 

Figure  10  shows  the  effects  produced  by  the  absorbing  material . 
The  dashed  curve  in  Figure  10  represents  the  E  plane  pattern  of  the 
apparatus  shown  in  Figure  9,  and  the  solid  curve  is  the  E  plane 
pattern  of  the  horn  and  ground  plane  alone. 

The  validity  of  this  rodded  medium  model  is  a  function  of  the 
specific  antenna  geometry  being  simulated.  If  the  Y  components 
of  the  induced  currents  are  producing  the  dominant  pattern  changes, 
the  rodded  medium  described  above  is  adequate  for  the  simulation. 

In  principle,  when  the  horn  aperture  is  large  compared  with  a 
wavelength  and  the  thickness  of  the  sheath  is  small  compared  to 
the  aperture,  the  majority  of  the  oscillatory  currents  in  the  plasma 
can  be  assumed  to  be  directed  along  the  Y  axis.  These  arguments 
are  valid  when  the  angle  6  is  not  too  large.  The  Z  components 
of  the  currents  become  increasingly  important  as  6  approaches 
90  degrees  in  the  E  plane. 
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simplified  side  view  of  rodded 

MEDIUM  AND  HORN  GEOMETRY 
(  RODS  IN  Y  DIRECTION  ) 


FIGURE  7  -  Antenna  geometry  with  rods  oriented  in  the 


direction. 
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Figure  8.  Experimental  antenna  geometry  without  the  absorbing 
material  (rods  in  the  Y  direction). 
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Figure  9 


Experimental  antenna  configuration  with  the 
absorbing  material  (rods  in  the  Y  direction). 
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NO  LOAD  PATTERNS  AT  II  KMC  SHOWING  THE 


EFFECTS  OF  ADDING  THE  ABSORBER  MATERIAL. 


-  E  PLANE  PATTERN 

WITHOUT  ABSORBER 


. - E  PLANE  PATTERN 

WITH  ABSORBER 


„  FIGURE  10  -  Radiation  patterns  with  and  without  the  absorbing 
■aterial. 
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Since  taere  are  currents  flowing  in  the  ^direction  in  this 

antenna  geometry,  an  improvement  in  the  simulation  of  a  homogeneous 

isotropic  plasma  layer  was  accomplished  by  orienting  the  wires 

in  both  the  Y  and  2  directions.  The  wires  parallel  to  the  X  axis 

were  omitted  on  the  premise  that  the  currents  excited  in  the 

X  direction  prod' ce  little  or  no  effects  in  the  F,  *nd  li  plane 

measurements.  The  wire  mesh  structures  used  in  the  fabrication 

of  this  artificial  dielectric  were  held  in  position  by  styrofoam 

spacers  (see  Figure  11).  Similar  techniques  have  been  utilized 
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in  the  construction  of  artificial  dielectric  lenses  .  The 
thickness  of  this  slab  was  2.7  inches  ,and  the  equivalent  plasma 
frequency  was  8.1  KMC.  The  experimental  antenna  is  shown  in  Figure  12. 

Pattern  Measurements 

The  radiation  patterns  were  measured  on  a  laboratory  apparatus 
shown' in  Figure  13.  The  probe  antenna,  which  was  supported  by  a 
movable  plywood  arm,  was  operated  as  a  receiving  antenna.  The 
transmission  path  between  the  transmitting  and  receiving  antennas 
was  6  feet  and  the  angular  position  of  the  movable  arm  with  respect 
to  the  axis  of  the  transmitting  horn  could  be  varied  from  0  degrees 
to  90  degrees  in  5  degree  increments.  It  became  apparent  that 
the  angular  resolution  of  this  apparatus  was  inadequate  for  this 
particular  experiment  and  some  of  the  fine  structure  of  the  patterns 
was  not  detected. 

In  order  to  obtain  better  angular  resolution,  the  measurements 
were  performed  on  a  free  space  antenn.  range  shown  in  Figure  14. 

The  antenna  positioner  was  equipped  to  rotate  in  both  the  9  and  $ 
directions.  The  transmitting  source  or  probe  horn  was  mounted  on 
a  stationary  non-reflecting  pole  located  about  30  feet  from  the 
receiving  antenna  at  a  height  of  17  feet.  The  free  space  range  was 
instrumented  to  measure  the  receiving  patterns  of  an  antenna  system. 
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THE  ARTIFICAL  DIELECTRIC  SHOWN  IN  THIS  DIAGRAM 
IS  CONSTRUCTED  WITH  THE  WIRE  MESH  STRUCTURE 
SANDWICHED  BETWEEN  DIELECTRIC  FOAM  SPACERS. 

THE  FOAM  SPACERS  ARE  GLUED  TOGETHER,  THEREBY 
HOLDING  THE  WIRES  IN  PLACE.  IT  SHOULD  BE  NOTED 
THAT  THE  WIRES  NEED  NOT  MAKE  ELECTRICAL 
CONTACT  AT  JUNCTIOF  POINTS 


FIGURE  11  -  Rodded  medium  with  rods  in  both  the  Y  and  2  directions. 
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Figure  12 . 


Experimental  antenna  geometry  using  the  rodded  medium 
shown  in  Figure  11  (rods  in  both  the  Y  and  Z  direction). 
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FIGURE  13  -  Laboratory  pattern  range 


ANT  SYSTEM  BEING  TESTED 
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FIGURE  14  -  Block  diagram  of  the  free  space  antenna  range 


However,  if  the  probe  horn  and  antenna  being  tested  are  loosely 

coupled  and  the  driving  point  impedances  are  independent  of  angular 

position,  by  reciprocity  the  receiving  and  transmitting  patterns 
28 

are  equal.  Figure  15  illustrates  the  comparison  between  the 
patterns  measured  with  the  laboratory  apparatus  and  a  free  space 
antenna  range. 


The  E  and  H  Plane  patterns  of  the  horn  antenna  with  and  without 
the  artificial  dielectric  were  measured  at  9,  10  and  11  KMC.  Both 
rodded  media  were  used  in  the  measurements.  The  patterns  obtained 
with  the  apparatus  shown  in  Figure  9  are  displayed  in  Figures  16 
through  22.  The  two  slab  thicknesses  chosen  for  these  tests  were 
1.5  inches  and  2.1  inches.  The  experimental  measurements  were 
compared  with  the  theoretical  calculations  of  the  appropriate  slit 
problem  referred  to  in  the  introduction.  The  E  piane  patterns  were 
compared  with  Equation  (1)  and  the  H  piane  patterns  were  compared 
with  Equation  (2). 

Equations  (1)  and  (2)  can  be  factored  into  two  functions,  one 
being  associated  with  the  pattern  without  the  plasma  sheath  and 
the  other  being  associated  entirely  with  the  plasma  sheath.  The 
resulting  patterns  can  be  expressed  as 
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4  =  thickness  of  the  sheath.  (S3) 

The  functions  and  L0  are  respectively  the  t  plane  and  H  plane 

patterns  of  a  horn  in  an  infinite  ground  plane.  Figures  23  and 

24  summarize  the  experimentally  determined  values  of  A  A*  and 

9  r 

the  corresponding  theoretical  computation  obtained  from  Equation  (50). 

* 

The  experimental  values  of  A  A  are  obtained  by  dividing  the 

9  9 

composite  H  plane  pattern  of  the  horn  and  artificial  dielectric  by 
the  H  plane  pattern  of  the  horn  alone. 

The  E  plane  patterns  obtained  with  the  apparatus  shown  in 
Figure  12  are  given  in  Figures  25  and  26,  The  artificial  dielectric 
used  in  these  measurements  had  rods  oriented  in  both  the  Y  and  Z 
directions,  and  the  equivalent  thickness  of  the  simulated  plasma  sheath 
was  2.7  inches.  Again  the  experimental  measurements  were  compared 
with  Equation  (1). 


i 
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NOTES-  (Of*  11.3 KMC 
8.1  KMC 

(3)X  *  1.5  inches 

_ FREE  SPACE  RANGE 

OOO  LABORATORY  RANGE 


FIGURE  15  -  Comparison  between  the  radiation  patterns  measured 
*  with  the  laboratory  and  free  space  pattern  ranges. 
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jt  -  0 ,  no  plosmo,  f=9KMC 

— EXPTL  DATA 
- THEORY  (REF  2,n«,426) 


FIGURE  16  -  H  plane  radiation  patterns  of  the  antenna  shown  in 
Fig.  9  at  9  KMC  (rods  in  the  Y  direction). 
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JL  =0  ,  f=9KMC,  NO  PLASMA 

-  EXPTL.  DATA 


- THEORY (REF  2, n  =  426) 


FIGURE  17  -  H  plane  radiation  patterns  of  the  antenna  shown  in 
Fig.  9  at  9  KMC  (rods  in  the  Y  direction). 


-£  =  l.5",f  =IOKMC,  fp=  8. IKMC 


sX=  0,  f  =  10  KMC,  NO  PLASMA 

- EXPTL. DATA 

- THEORY(REF.  2tn  =  564) 


FIGURE  18  -  H  plane  radiation  patterns  of  the  antenna  shown 
in  Fig.  9  at  10  KMC  (rods  in  the  Y  direction). 
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- EXPTL.  DATA 

- THEORY  (REF  2,  n  =  673) 


FIGURE  19  -  H  plane  radiation  patterns  of  the  antenna  shown  in 
Fig.  9  at  11  KMC  (rods  in  the  Y  direction). 


J,  s  0  %  NO  PLASMA,  f  =  9  KMC 

- EXPTL.  DATA 

- THEORY  (  REF  2,  D  =  .426) 

FIGURE  20  -  E  plane  radiation  patterns  of  the  antenna  shown  in 
Fig.  9  at  9  KMC  (rods  in  the  Y  direction). 
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J?=  0"  NO  PLASMA,  f  =  10  KMC 

- EXPTL.  DATA 

- THEORY (  REF  2,  n=  .564) 


FIGURE  21  -  E  plane  radiation  patterns  of  the  antenna  shown  in 
Fig  9  at  10  KMC  (rods  in  the  Y  direction). 


i=0,  NO  PLASMA,  f  =  II  KMC 


- EXPTL.  DATA 

- THEORY  (REF  2, n«.675) 


FIGURE  22  -  E  plane  radiation  patterns  of  the  antenna  shown  in 
Fig.  9  at  11  KMC  (rods  in  the  Y  direction}. 
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PLANE  DATA  AT  9 KMC  (INDEX  OF  REFRACTION  =  426  AND  THE 

EQUIVALENT  PLASMA  FREQUENCY  =  8.1  KMC) 


FIGURE  23  -  Measured  values  of  A  ,  A  ,  at  9  KMC. 


FIGURE  24  -  Measured  values  of  A 


-&=0,  NO  PLASMA,  f  =10  KMC 

- EXPTL.  DATA 

- - THEORY  (REF  2,  H  =  .564) 


FIGURE  25  -  E  plane  radiation  patterns  of  the  antenna  shown  in 

Fig.  12  at  10  KMC  (rods  in  both  the  Y  and  Z  directions) 


j£=0,  NO  PLASMA,  f  =11  KMC 


-  EXPTL.  DATA 

- THEORY  (  REF  2 ,  D  =  .675  ) 

FIGURE  26  -  E  plane  radiation  patterns  of  the  antenna  shown  in 

Fig.  12  at  11  KMC  (rods  in  both  the  Y  and  Z  directions). 


Discussion  of  Experimental  Results 

The  correlation  between  the  H  plane  measurements  and  the 
theoretical  calculations  given  by  Equation  (2)  is  very  good,  as 
seen  from  Figures  23  and  24.  However,  the  agreement  between  the 
E  plane  measurements  and  the  theoretical  computations  is  less 
satisfactory.  The  sharp  peaks  in  the  radiation  pattern  occurring 
near  the  critical  angle 


sin  e„  =  n, 
c 

predicted  by  Equation  (1),  were  not  observed  in  the  E  plane  measure¬ 
ments.  These  sharp  resonance  phenomena  require  the  correct  phase 
relationships  between  all  the  radiating  sources,  and  any  deviation 
from  this  coherent  phase  relationship  will  tend  to  destroy  the 
resonance.  The  lack  of  agreement  between  theory  and  experiment 
is  partly  attributable  to  the  differences  between  the  theoretical 
and  experimental  geometries.  The  experimental  measurements  are 
performed  with  a  horn  antenna  and  a  simulated  plasma  slab  of  finite 
size,  while  the  theoretical  model  is  based  on  an  infinite  slot 
antenna  and  an  unbounded  plasma  layer.  The  rodded  medium  also 
exhibits  anisotropic  behavior,  and  relative  rod  spacings  d/XQ 
used  in  these  experiments  are  approaching  an  upper  limit  where 
the  artificial  dielectric  starts  to  deviate  from  the  behavior  of 
a  simple  dielectric  material.  Aside  from  the  lack  of  agreement 

around  the  critical  angle  0  ,  the  radiation  patterns  are  predicted 

2C 

reasonably  well  by  Omura's  analysis. 

The  addition  of  the  wires  in  the  Z  direction  seems  to  improve 
the  agreement  between  the  measured  and  computed  patterns  to  some 
extent,  but  the  very  sharp  resonance  is  still  not  observed  (see 
Figures  21,  22,  25  and  26.) 
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Since  the  physical  geometries  of  the  two  rodded  media  arc  some¬ 
what  different,  a  quantitative  statement  about  the  improvements 
observed  by  adding  the  set  of  wires  perpendicular  to  the  ground 
plane  is  not  warranted.  The  antenna  geometry  used  with  the  rodded 
medium  having  rods  in  the  Y  direction  only  is  shown  in  Figure  9,  and 
the  antenna  geometry  used  with  the  artificial  dielectric  having  rods 
in  both  the  Y  and  Z  directions  is  shown  in  Figure  12.  In  one 
case  there  were  plastic  rods  and  absorbing  material  present,  and  in  the 
other  case  only  the  horn,  ground  plane,  and  rodded  medium  were  present. 
In  addition,  one  of  the  simulated  slabs  had  a  thickness  of  1.5  inches 
while  the  other  had  a  thickness  of  2.7  inches.  Further  experimental 
investigations  are  needed  in  this  area  to  gain  a  better  understanding 
of  the  behavior  of  the  different  rodded  media  configurations. 

In  general,  the  radiation  from  the  horn  and  simulated 'plasma 
sheath  exhibits  a  rapid  attenuation  at  angles  greater  than  the 
critical  angle  Gc,  snd  the  functions  Aq  Aq  and  A<j,  A'J  have  a 
maximum  occuring  at  an  angle  slightly  smaller  than  0C.  The  net 
result  is  that  the  radiation  patterns  may  be  drastically  modified 
by  the  presence  of  a  plasma  layer,  as  illustrated  in  Figures  16  and  17. 

3  4 

Flock  *  predicted  that  the  radiation  pattern  would  have  a 
secondary  maximum  due  to  the  radiating  sources  in  the  plasma,  but 
the  position  of  the  maximum  is  in  disagreement  with  the  experi¬ 
mental  measurements  because  of  the  assumptions  used  in  the 
theoretical  model.  Better  correlation  would  be  expected  if  more 
realistic  current  distributions  in  the  plasma  layer  were  employed 
in  the  analysis. 
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IV. 


CONCLUSIONS 


This  study  has  shown  that  the  tensor  permittivity  of  a  homo¬ 
geneous  plasma  in  an  infinite  D.C.  magnetic  field  can  be  simulated 
by  a  rodded  medium.  The  rods  in  the  artificial  dielectric  simulate 
the  effects  produced  by  the  motions  of  the  electrons  in  the  plasma. 
Under  the  influence  of  an  infinite  magnetostatic  field,  the 
electrons  in  the  plasma  are  constrained  to  move  along  the  magnetic 
lines  of  force;  similarly,  the  currents  induced  in  the  rodded  medium 
arc  constrained  to  flow  in  the  metallic  wires  thereby  simulating  the 
effects  produced  by  the  free  electrons  in  the  plasma.  Waveguide 
measurements  confirmed  the  anisotropic  behavior  of  the  rodded 
medium  when  the  rods  are  oriented  in  one  direction  only. 

The  pattern  experiments  further  substantiate  the  isotropic 
behavior  of  the  rodded  medium  in  a  free  space  geometry.  The  antenna 
pattern  measurements  agree  reasonably  well  with  the  theoretical 
model  given  in  Reference  2,  and  demonstrate  the  effects  a  plasma 
sheath  can  exert  on  an  antenna  system.  The  pattern  distortion 
observed  in  this  experiment  illustrates  another  important  character¬ 
istic  of  plasma  sheath  transmission  in  addition  to  the  more  fully 
studied  attenuation  effects. 

With  reference  to  the  horn  geometry,  it  was  shown  that  the 

simulation  of  the  homogeneous  isotropic  properties  of  a  plasma 

sheath  could  be  achieved  by  a  rodded  medium  which  has  wires  oriented 

in  both  the  Y and  Z directions.  It  was  further  observed  that  the 

omission  of  wires  in  the  Z direction  altered  the  L  plane  pattern 

in  the  neighborhood  of  the  critical  angle  0  and  had  little  effect 

c 

on  the  ti  plane  patterns. 
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V.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


s 


More  experimental  studies  on  the  anisotropic  behavior  of  the 
different  rodded  media  configurations  are  needed.  It  has  been  shown, 
from  a  theoretical  standpoint,  that  a  rodded  medium  could  simulate 
the  anisotropic  properties  of  a  cold  plasma  in  an  infinite  magneto¬ 
static  field,  but  a  definitive  experimental  study  would  indicate  to 
what  degree  the  simulation  is  taking  place.  Experimental  studies 

in  the  areas  of  eliminating  re-entry  blackout  by  applying  a  large 

29  30  31 

U.C.  magnetic  field  *  *  ,  the  magnetic  window  effect,  could  be 

demonstrated  using  an  anisotropic  rodded  medium  with  the  rods 
oriented  in  one  direction. 


The  addition  of  a  loss  mechanism  to  the  artificial  dielectric 
should  also  be  investigated  with  the  objective  of  developing 
techniques  to  simulate  both  the  plasma  and  collisional  frequencies  in 
a  free  space  geometry.  The  lossy  artificial  dielectric  could  be  used 
to  simulate  more  realistic  plasma  sheaths  in  laboratory  antenna 
studies. 


Collisional  effects  can  be  added  to  the  rodded  medium  by 
employing  lossy  wires.  Rotman^  showed  that  the  propagat  •m  constant 
of  a  rodded  medium  using  lossy  rods  could  be  expressed  as 

Z 

cosh  yd  =  cos  k  d  ♦  i  sin  k  d 

o  at  0 

where 


y  =  a  ♦  jB  =  complex  propagation  constant 

of  the  rodded  medium 


The  electric  field  is  assumed  to  be  parallel  with  the  rods  and  Rs 
is  the  internal  surface  resistance  of  one  of  the  equivalent  wire 
grating  structures.  The  surface  resistance  is  equal  to 
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where  RQ  is  the  resistance  per  unit  length  of  one  of  the  rods. 

Ni-chrorae  wire  can  be  utilized  in  the  construction  of  a  lossy  rodded 

medium,  but  the  equivalent  collisional  frequency  is  small.  Larger 

and  more  realistic  collision  frequencies  can  be  obtained  by  using  rods 

which  are  constructed  out  of  glass  or  dielectric  fibers  and  ar~ 

32 

coated  with  a  thin  metallic  film  .  The  resistance  per  unit  length, 
Rq,  of  the  dielectric  coated  rods  can  be  of  the  order  of  hundreds  of 
ohms  per  inch  instead  of  the  low  values  obtainable  with  the  Ni-chrome 
wire. 


Plasma  simulation  could  also  prove  to  be  useful  in  microwave 
diagnostic  studies  of  laboratory  plasmas.  The  determination  of 
some  of  the  experimental  limitations  of  a  microwave  diagnostic 
geometry  (such  as  the  errors  introduced  by  the  finite  experimental 
geometry)  could  be  obtained  by  replacing  the  actual  laboratory 
plasma  by  a  rodded  medium. 

A  more  practical  application  of  the  anisotropic  characteristics 

of  a  rodded  medium  is  the  construction  of  a  transmission  type 

polarizer  which  converts  a  plane  polarized  wave  into  an  ellipticallv 
33 

polarized  wave  .  Tf  the  rods  are  oriented  at  45°  to  the  incident 
electric  field  of  a  plane  wave  and  the  parameters  of  the  slab  are 
chosen  correctly,  the  emerging  wave  will  be  circularly  polarized. 
This  type  of  wave  filter  in  optics  is  commonly  called  a  "quarter- 
wave  plate."  'Hie  development  of  a  quarter-wave  plate  which  would 
convert  a  linear  polarized  antenna  into  a  circularly  polarized 
antenna  could  prove  to  be  extremely  useful. 


-54- 


REFERENCES 


1. 


2. 


j, 


4. 


5. 


6. 


7. 


3. 


9. 


10. 


11. 


12. 


Tamir,  T.  and  A.  A.  Oliner,  "The  influence  of  complex  waves  on 
the  radiation  field  of  a  slot-excited  plasma  layer," 
IRE  Trans.  PGAP,  vol.  AP-10,  pp.  55-65,  January,  1962. 

Omura,  M. ,  "Radiation  pattern  of  a  slit  in  a  ground  plane 
covered  by  a  plasma  layer,"  Air  Force  Cambridge 
Research  Laboratory,  AFCRL-62-958,  December, T962. 

Flock,  W.  L,,  Plasma  effects  on  microwave  radiation,  Ph.D. 
dissertation,  Dept,  of  Engineering,  University  of 
California,  Los  Angeles,  June,  1960. 

Flock,  W.  L.  and  R.  S.  Elliott,  "The  radiation  pattern  of  a 
microwave  horn  and  a  plasma  layer,"  IRE  Trans.  PGAP, 
vol.  AP-10,  pp.  65-88,  January,  1962. 

Clayton,  L.  and  J.  S.  Hollis,  "Calculation  of  microwave 

antenna  radiation  systems  by  the  Fourier  integral 
method,"  Microwave  Journal ,  vol.  3,  no.  9,  pp.  59-66, 
September,  1960. 

Rotman,  W.,  "Plasma  simulation  by  artificial  dielectric  and 
parallel-plate  media,"  IRE  Trans.  PGAP,  vol.  AP-10, 
pp.  82-95,  January,  1962. 

Collin,  R.  E.,  "A  simple  artificial  anisotropic  dielectric 
medium,"  IRE  Trans.  PGMTT ,  vol.  MTT-6,  pp.  206-209, 
April,  1958. 

Kirschbaum,  H.  S.  and  S.  Chen,  "A  method  of  producing  broad¬ 
band  circular  polarization  employing  an  anisotropic 
dielectric,"  IRE  Trans.  PGMTT,  vol.  MTT-5,  pp.  199-203, 
July,  1957. 

Estrin,  G. ,  "The  effects  of  anisotropy  in  a  three  dimensional 
array  of  conducting  disks,"  Proc.  IRE,  vol.  39, 
pp.  821-826,  July,  1951. 

Delcroix,  J.  L. ,  Introduction  to  the  theory  of  ionized  gases. 
Interscience  Publishers,  New  York,  1960,  pp.  67-73. 

Elliott,  R.  S.,  "The  tensor  permittivity  of  a  weakly-ionized 
plasma,"  Microwave  Journal,  vol.  2,  no.  3, 
pp.  47-50,  Fiarch",  1959. 

Whitmer,  R.  F.,  "Principles  of  microwave  interactions  with 
ionized  media,"  Microwave  Journal,  vol.  2,  no.  3, 
pp.  47-50,  March,  1959. 


-5b- 


REFERENCES  (Continued) 

13.  Brown,  M.  A,,  "Artificial  dielectrics  having  refractive  indices 

less  than  unity,"  Proc.  IBS,  Monograph  no.  62R, 
volt  100,  part  4,  pp.  5 1-62 ,  May,  1953. 

14.  Brown,  J.  and  W.  Jackson,  "The  properties  of  artificial 

dielectrics  at  centimeter  wavelengths,"  Proc.  IBB, 
paper  no.  1699R,  vol.  102B,  pp.  11-21,  January,  1955. 

15.  Cohn,  S.  B.,  "Analysis  of  the  metal  strip  delay  structure  for 

microwave  lenses,"  J.  Appl.  Phys.,  vol.  20, 
pp,  257-262,  March,  Y<R9. 

16.  Cohn,  S.  B.,  "Blectrolytic-tat.  t  measurements  for  microwave 

metallic  delay  lens  media,"  J.  Appl.  Phys.,  vol.  2.1, 
pp.  674-680,  July,  1950. 

17.  Collin,  R.  E.,  Field  theory  of  guided  waves,  McGraw-Hill, 

New  York,  1960,  chapt.  9. 

18.  Brown,  J.,  "The  design  of  metallic  delay  dielectrics," 

Proc,  IEB,  vol.  97,  part  3,  pp.  45-48,  January,  1950. 

19.  Macfarlane,  G.  G. ,  "Surface  impedance  of  an  infinite  parallel- 

wire  grid  at  oblique  angles  of  incidence," 

Journal  IBB,  vol.  93,  part  3A,  pp.  1523-1527,  1947. 

20.  Marcuvitz,  N.,  Waveguide  handbook,  McGraw-Hill,  New  York, 

1951,  p.  285. 

21.  Jasik,  H. ,  Antenna  engineering  handbook,  McGraw-Hill. 

New  York,  1961,  chapt.  14,'  p.  23 . 

22.  Kittle,  C.,  Introduction  to  solid  state  physics,  Wiley, 

New  York ,  1956 ,  pp.  162-163.  *— — — 

23.  Von  Hippie,  A.,  Dielectrics  and  waves,  Wiley,  New  York, 

1954,  pp.  97-98. 

24.  Moon,  P.  and  L).  E.  Spencer,  Foundations  of  electrodynamics. 

Van  Nostrand  Co.,  Princeton,  N.  J.,  1960,  p.  199. 

25.  Montgomery,  C.  G. ,  Technique  of  microwave  measurements, 

McGraw-Hill,  New  York,  1947,  pp.  565-566. 

26.  Harrington,  R.  F. ,  Time-harmonic  electromagnetic  fields, 

McGraw-Hill,  1961,  p.  67. 

27.  Kock,  W.  E.,  "Metallic  delay  lenses,"  Bell  Systems  Technical 

Journal,  vol.  27,  pp.  58-82,  January,  1948. 


-56- 


REFERENCES  (Continued) 


28.  Silver,  S. ,  Microwave  antenna  theory  and  design,  McGraw-Hill, 

New  York,  I9T9,  pp.  48^49. 

29.  Dirsa,  E.  F.,  "The  telemetry  and  communication  problems  of 

re-entrant  space  vehicles,"  Proc.  IRE,  vol.  48, 
pp,  703-713,  April,  1960. 

30.  Hodara,  H. ,  "The  use  of  magnetic  fields  in  the  elimination  of 

the  re-entry  blackout,"  Proc,  IRE,  vol.  49, 
pp.  1825-1830,  December,  1961. 

31.  Hodara,  H.  and  G.  I.  Cohn,  "Radiation  characteristics  of  a 

slot  antenna  in  the  presence  of  a  lossy  anisotropic 
plasma  sheath,"  Proc,  of  the  National  Electronics 
Conference,  Chicago,  Ill.,  October  1961,  vol.  17, 
pp.  622-631 . 

32.  Golden,  K.  E.  and  T.  M.  Smith,  ‘Simulation  of  a  thin  plasma 

sheath  by  a  plane  of  wires,"  To  be  published  in 
IEEE  Transactions  on  Nucl.  Sci.,  vol.  NS-11,  no.  1, 
January,  1964. 

33.  Ploussios,  G. ,  "A  circular  polarizer  for  linear  polarized 

antennas,"  Paper  presented  at  the  Northeast 
Electronics  Research  and  Engineering  Meeting,  Boston, 
Mass.,  November  4-6,  1963. 

34.  Jackson,  J.  D.,  Classical  electrodynamics,  Wiley,  New  York, 

1962,  p.T 

35.  Stone,  J.  M. ,  Radiation  and  optics,  McGraw-Hill,  New  York, 

1963,  cHaptT  17. 

36.  Koerbcr,  G.  G.,  Properties  of  solids,  Prentice-Ilall, 

Englewood  Cliffs,  N.  J.,  1962,  chapt.  7. 

37.  Stratton,  J.  A.,  Electromagnetic  theory,  McGraw-Hill, 

New  York,  1941,  p.  66. 


-57- 


APPENDIX  A 


Plane  Wave  Propagation  in  Anisotropic  Dielectric 

34  35 

Plane  wave  type  of  propagation  in  an  anisotropic  dielectric,  * 
will  be  discussed  in  this  section.  The  propagation  characteristics 
for  a  monochromatic  wave  with  a  time  dependence  given  by  e+^ut 
are  obtained  from  Maxwell's  Equations  . 


V  X  E  =  -  jup  U 

da) 

V  X  H  =  ju  U 

(2a) 

U  =  e0  [er]  [E] 

(3a) 

Assuming  the  tensor  | e r |  to  be  diagonalized  with  the  coordinate  axes 
oriented  along  the  principal  directions  of  the  dielectric,  Equation 
(3a)  becomes 


U1  =  eo  ell  L1 

U2  =  Co  C22  L2 

D3  =  eo  e33  L3  * 

Taking  the  curl  of  Equation  (la) 

V  x  V  x  E  =  -  juip  V  x  II 

J  o 

unu  replacing  V  X  H  with  jwb  yields  the  equation  governing  the 
propagation  characteristics  of  the  medium. 

V  x  V  x  E  -  k~lcJ[E]=0 

O  T  1  1 


(4a) 

(5a) 

(6a) 


(7a) 


(8a) 
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-  i  k  *  r 

If  the  propagation  of  the  form  e  is  assumed,  the  curl  of  the 
curl  of  the  electric  field  can  be  replaced  by 

V  x  7  x  I;  *  -  k  x  k  X  L  .  (9a) 

Let  the  vector  wavenumber  k  be  renresented  by 

k  =  kn  ,  (10a) 

where  n  is  a  unit  vector  in  the  direction  of  k.  The  substitution 
of  Equations  (9a)  and  (10a)  into  Equation  (8a)  yields  the  following 
expression: 


n  x  n  x  E  ♦  k^/k2  (er)  [E]  =  0 

The  vector  identity 

n  x  n  x  E  =  (n«E)  n  -  E 

is  used  to  simplify  Equation  (11a). 

(n«E)  n  -  E  ♦  k2/k2  [er]  [E]  =  0 

Ihe  rectangular  components  of  Equation  (13a)  are 


(11a) 


(12a) 


(13a) 


.2,  .2 
kQ/k 

£n  * 

■!•>] 

|  E1  ♦  n2  nl  h  +  "3  nl 

U3  »  0 

(14a) 

nl  n2 

Li  * 

2  .  | 
e22  ♦  n2  -  !j 

C2  *  n3 

n2  i;3  =  0 

(15a) 

“l  n3 

hi  * 

n2  n3  1 

t.2.,2  2 

r_  ♦  1  k  /k  e,,  ♦  n,  -  1 
2  1  o  33  3 

Es  '  0 

(16a) 

where  k  is  the  free  space  wavenumber  and  n, ,  n^,  n~  are  respectively 
tiie  Xj,  and  X.,  components  of  the  unit  vector  n.  The  solutions 
for  k^  are  obtained  by  setting  the  determinant  of  the  above  set  of  • 
homogeneous  equations  equal  to  zero. 
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In  general,  there  are  two  modes  of  propagation  described  by 

Equations  (14a),  (15a)  and  (16a).  One  mode  has  a  Poynting  vector 

which  is  not  parallel  with  the  direction  of  the  wave  motion,  and 

the  other  mode  has  a  Poynting  vector  and  vector  wave  number  in  the 

2 

same  direction.  Table  (la)  summarizes  the  values  of  k  for  different 
polarizations  and  different  directions  of  propagation. 

TABLE  la 


Polarization  of  Direction  of  Propagation  Constant 

the  Electric  Field  Propagation _ 

X.  Axis  X-  X_  Plane  k2  =  k2  c,. 

1  4  3  O  11 


2  2 

X~  Axis  X.  X,  Plane  k  » ■■  k 

2  13  O  22 

X,  Axis  X.  X_  Plane  k2  -  k2  c. 

3  12  •  o  j3 

Xj  X2  Plane  Xx  X£  Plane  k2  =  k2 

ell  C22 

T"  2 

\  hi  *  n2  e22 

X.  X.  Plane  X.  X.  Plane  k2  =  k2 

13  13  0 

'll  e33 

*T“.  “  2 

-nl  hi  n3  e33 

X-  Plane  X0  X_  Plane  k2  -  k2 

2  3  2  3  o 

e22  e33 

2  2 
_n2  e22  *  n3  e33 
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APPENDIX  B 


Certain  restrictions  are  placed  on  the  components  of  the  tensor 
permittivity  of  the  rodded  medium  because  of  the  symmetries  of  the 
structure.  In  general,  the  tensor  permittivity  of  an  anisotropic 
material  is  expressed  by 


D. 


D. 


C11  e12  e13 


e21  e22  e23 


e31  e32  C33 


E, 


E. 


E. 


(lb) 


36 

The  constraints  arc  obtained  by  applying  Neumann's  principle 
when  a  symmetry  operation  is  performed.  In  other  words,  the  components 
remain  invariant  under  a  coordinate  transformation  which  corresponds 
to  a  symmetry  operation. 


The  rodded  mediuta  under  discussion  in  this  section  has  all  the 
rods  oriented  in  the  X^  direction  forming  an  infinite  array  of  con¬ 
ducting  cylinders  throughout  all  space.  The  following  symmetry 
operations  can  be  imposed  on  the  rodded  medium  because  of  the 
structural  symmetries  shown  in  Figure  2. 

(1)  180  degree  rotation  about  the  axis. 

(2)  180  degree  rotation  about  the  X2  axis. 

(3)  90  degree  rotation  about  the  X^  axis. 


Let  the  unprimed  quantities  be  referred  to  the  original  coordinate 
system  and  the  primed  quantities  be  referred  to  the  transformed 
coordinate  system.  The  coordinate  transformation  can  be  expressed  by 
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11 

a12 

a13 

21 

a22 

a23 

31 

rsi 

rO 

c4 

a33 

(2b) 


The  components  of  [e]  in  the  primed  coordinate  system  are  given  by' 


3  3 

l  l 

i*l  j-1 


a  ■  a  .  e.  . 
mi  nj  ij 


(3b) 


where  a  .  and  a  .  are  the  components  of  the  transformation  matrix 
mi  nj  1 

[A]  shown  in  Equation  (2b) .  The  transformation  matrix  for  a  180  degree 
rotation  about  the  Xj  axis  is  equal  to 


10  0 
0-10 
0  0-1 


(4b) 


The  components  of  the  tensor  permittivity  obtained  from  Equation  (3b) 
are  listed  below. 


11 

** 

eii 
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12 
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"e12 
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13 

s 

”e13 

i 

21 

= 

“e21 

i 

22 

3 

e22 

i 

23 

3 

e23 

i 

31 

M 

“e31 

i 

32 

= 

e32 

i 

33 

m 

£33 
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The  application  of  Neumann's  principle 


» 


c 


ij 


yields  the  following  allowable  components  of  the  tensor  permittivity. 


C11 

0 

0 

Id  -  e0 

0 

e22 

e23 

(Sb) 

0 

e32 

c33 

The  transformation  matrix  for  the  180  degree  rotation  about  is 
given  by 


[A] 


-10  0 
0  10 
0  0-1 


(6b) 


and  allowable  components  of  [c]  under  this  symmetry  operation  are 


M 


en 

o 

e31 


0 

e22  0 

0  e 


13 


33 


(7b) 


The  third  symmetry  operation ,  a  90  degree  rotation  about  the 
axis,  is  given  by  the  transformation  matrix 


IA] 


0  1  0 

-10  0 
0  0  1 


(8b) 


-65- 


and  imposes  the  following  constraints  on  the  values  of  c^. 


‘■a 

e12 

e13 

[e]  * 

'  e12 

ell 

e13 

e31 

e31 

e33 

(9  b) 


Since  Equations  (5b) (  (7b)  and  (9b)  must  be  satisfied  simultaneously, 
the  following  components  are  zero. 


c 

e 

e 

c 

c 

e 


12 

13 

21 

23 

31 

32 


=  C 

«  0 

a  0 
«  0 
*  0 
a  0 


Therefore,  the  resulting  tensor  permittivity  of  the  rodded  medium 
illustrated  in  Figure  2  is  equal  to 


M 


ell 

0 
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0 

E11 

0 


0 

0 


(10b) 
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